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, ABSTRACT 

Atmospheric temperature and mixing ratio profiles of terrestrial planets vary with the spectral en- 
f — ■ ergy flux distribution for different types of M-dwarf stars and the planetary gravity. We investigate 

the resulting effects on the spectral appearance of molecular absorption bands, which are relevant 
as indicators for potential planetary habitability during primary and secondary eclipse for tran- 
£\j | siting terrestrial planets with Earth-like biomass emissions. Atmospheric profiles are computed 

using a plane-parallel, ID climate model coupled with a chemistry model. We then calculate 
simulated spectra using a line-by-line radiative transfer model. 

We find that emission spectra during secondary eclipse show increasing absorption of 
methane, water, and ozone for planets orbiting quiet M0-M3 dwarfs and the active M-type star 



o 



X 



■ AD Leo compared with solar-type central stars. However, for planets orbiting very cool and quiet 

M dwarfs (M4 to M7), increasing temperatures in the mid-atmosphere lead to reduced absorption 
signals, which impedes the detection of molecules in these scenarios. Transmission spectra dur- 
ing primary eclipse show strong absorption features of CH4, N2O and H2O for planets orbiting 
quiet M0-M7 stars and AD Leo. The N2O absorption of an Earth-sized planet orbiting a quiet 
M7 star can even be as strong as the CO2 signal. However, ozone absorption decreases for plan- 
ets orbiting these cool central stars owing to chemical effects in the atmosphere. To investigate 
the effect on the spectroscopic detection of absorption bands with potential future satellite mis- 
sions, we compute signal-to-noise-ratios (SNR) for a James Webb Space Telescope (JWST)-like 
aperture telescope. 

Key words. Planets and satellites: atmospheres; Astrobiology 
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1. Introduction 



The majority of the more than 500 extra-solar planets found so far are giant gas planets. However, 
with improving detection techniques, more and more low-mass planets will be found. At the time 
of writing this paper, we know of more than 20 planets with (minimum) masses below 1 M m . Two 



systems (around Gliese 581 and HP 40307) ho st three potential low-mass planets each (lUdrv et al 



2007 



Mavor et al 



2009b 



Mavor et al 



2009al). The only two low-mass planets wher e both radius 



and mass could be measured are CoRoT-7b (ILeger et al 



2009 



Ouel oz et al 



2009), a transiting 



2009), with 6.5 M e 



planet of approximately 5 M ffi and 1.7 R @ , and GJ 1214b (Charbon neau et al. 
and 2.7 R B . CoRoT-7b, which is likely to be a rocky planet, may therefore be classified as terrestrial 
(a so called super-Earth). 

Whereas atmospheres of close-in gas giant planets around bright stars have been already 



telescope (e.g. 


Vidal-Madiar et al.ho04 


, Tinetti et al. 


2007 


Swain et alj|2008 


Knutson et al.lho08 


Swain et al. 


2009 


Knutson et al. 


2009) 


, terrestrial extra-solar planets seem out of reach with cur- 



rently available telescopes. However, future satellite missions, such as the JWST (James Webb 
Space Telescope) or proposed dedicated missions for spectral characterisation, will provide an im- 
proved detection range. To which extent their performance will be sufficient to spectroscopically 
characterise the atmosphere of terrestrial planets in the habitable zone of their central star is the 
subject of increasing modelling investigations. 

Previous major modelling efforts pro ducing synthe ti c spectra o f hypo t hetical terrestrial 
extra- s o lar planets w er e perf or med by e.g 



Selsis 



(120021). ISeguraet al 



(2003), 



Segura et al 



(2000) 



(2005), 



Selsis et al 



Tinetti 



(2006), 



(2002) 



Des Marais et al 



Ehrenreich et al 



(2006), 



Kaltenegger et al 



(2007). 



Kaltenegger & Traubl (120091) and lKaltenegger & Sasselov 



20101) . They 



examined the influence of e.g. stellar type, atmospheric abundances, background atmospheres or 
atmospheric evolution on the emission and transmission spectra. 

This work focuses on emission and transmission spectra of super-Earth planetary atmospheres 
that orbit in the habitable zone (HZ) around M-dwarf stars because in these cases the highest signal- 
to-noise-ratios (SNRs) for potentially habitable planets are expected. Spectral resolution and the 
expected SNRs are studied for a photon-limited telescope with 6.5 m aperture (JWST-sized). Real 
SNRs will be significantly lower because of instrumental noise and the astronomical background. 
The values here serve as a "best case" therefore and provide a first estimate on the possible detection 
performance. A parametric study of expected SNRs for a wide ra nge of star-planet combinations 



Belu etal 



d2010l) . We comple- 



including instrumental noise for the JWST has been performed by 
ment this study by a detailed investigation that includes the effects of atmospheric chemistry on 
habitable super-Earth planets. 

Our approach is similar to that of 



Selsis 



(2000) 



Segura et al 



(2003J) and 



Segura et al 



(2005) 



We first apply an atmosphere column model to calculate the mean temperature profile and the cor- 
responding chemical profiles assuming an Earth-like planet development and biomass emissions. 
We chose this approach because the T-p-profiles as well as the atmospheric mixing ratios vary de- 
pending onstejj^rjype andplanetary parameters. This was first shown for planets orbiting M-dwarf 
stars bv lSegura et al.l (120051) . who demonstrated the effect of the changing atmospheric chemistry 
on emission flux spectra of these Earth-like extra-solar planets. This early work did not include 
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transit transmission spectra and SNR estimates or a cha nge of planetary gravi t y thou gh. SNRs were 
later calculated for a 6.5m JWST-aperture telescope bv lKaltenegger & Tra ub (2009) in a simplified 
approach, which kept the planetary atmospheric conditions fixed to those of a modern Earth when 
changing the central star. Their results showed the difficulty to detect the weak biosignatures in 
transmission spectra with the JWST, but also discussed how this could be overcome by co-adding 
several transit measurements. Our main goal is to extend these SNR investigations by including 
the effect of chemistry to study whether the detecti on of biosignatures could be favoured by chem- 
ical processes, as suggested bv lSegura et al.l (120051) . In addition, we also consider transmission and 
emission spectra. 

The paper is organised as follows: Section 2 describes the atmospheric and radiative transfer 
models and the modelled scenarios. In Sect. 3 we show how the SNR calculations are performed. 
Section 4 compares the resulting emission and transmission spectra with modern Earth spectra, 
followed by the resolution and SNR calculations. Section 5 summarises the results. 



2. Models and simulated scenarios 



2.1. Atmospheric and spectral models 

We used a ID radiative-convective photochemical model to calculate the chemistry and climate 
of planets with an Earth-like atmospheric chemical composition. The model calculates the glob- 
ally, diurnally-averaged atmospheric temperature, pressure, a nd composition profil es for cloud-free 
conditions. T he original code has been described in detail b y iKasting et al. ( 1984) and was further 



developed by ISegura et al 



(2003) 



Grenfell et al 



(l2007al) and 



Grenfell et al 



d2007bl) . Here, we 



briefly summarise the model and describe its recent improvements. 

The model extends from the ground up to the mid-mes osphere. The climate module grid is 



Grenfell et al 



d2007al) . we improved 



divided into 52 layers on a log-pressure grid. Compared with 
the calculation of the grid. Depending on how much water is evaporating at the surface, the pressure 
is re-adjusted in each iteration step. The chemistry module is divided into 64 equidistant layers 
in altitude. As an upgrade to previous model versions, the top layer is now equal to the top-of- 
atmosphere (TOA) altitude derived by the climate module. Hence, the model grid consistently takes 
into account the changing scale height with varying gravity in the climate and chemical modules, 
allowing the simulation of super-Earth planets with high gravity with consistent resolution. 

The radiative transfer for incoming stellar radiation in our climate module covers the wave- 
length region from 237.6 nm to 4.545 nm with 38 wavelength bands. We used a two-stream radia- 



tive transfer method based on 



Toon et al 



( 1989). The radiatively active species considered are H2O, 



CO2, O2, O3, and CH4. Absorption coefficients are calculated by an exponential sum correlated-k 
method using up to four terms per sum. Rayleigh scattering is included for N2, O2, and CO2. 

For thermal emission we considered the wavel ength range from 3. 33 - 1000 ^m in 16 bands, 



Mlawer et al.ll 1997 ). The model includes H2O, 



using the Rapid Radiative Transfer Model (RRTM, 
CO2, O3, CH4, and N2O as radiative gases in the IR-range. The absorption coefficients were cal- 
culated with the 16-term correlated-k method. Because k-coefficients are always calculated for a 
certain composition, temperature, and pressure range, the model is valid strictly only within 0.01 - 
1050 mbar and for the Earth's atmospheric temperature profile ±30 K. For this reason, we limit our 
coolest input stellar spectra up to M7 dwarf stars for Earth-sized planets, and up to M5-type stars 
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for super-Earth planets. Cooler M dwarfs would result in T-p profiles widely outside the applicable 
model range. Wet or dry adiabatic H2O lapse rates are used to o btain the temperature profile i n the 
troposphere where applicable. The water profile is taken from iManabe & Wetheraldl (119671) . We 
consider a mean solar zenith angle of 60° for the climate model. 

The chemistry module includes 55 species with 212 chemical reactions. As in most previous 
model approaches and for modern Earth simulations, CO2, O2 and Nt are isoprofiles. To calculate 
the photolysis rates, we took into account a wavelength range from 121.4 - 855 nm and considered 
O3, O2, CO2, H2O, SO2, CH4, and H2S as radiative gases. Rayleigh scattering was included for 



N2, O 2, and CO2. The mean solar zenith angle is 45°, following t he approach by ISegura et al 



(2005). In addition to the species considered by ISegura et al 



(12005b , we included H 2 and CH 4 



when calculating the mean molecular weight of air and the concentration of the fill gas nitrogen. 
We furthermore changed the convergence criterion for the numerical solution from previous model 
versions to a maximum difference in volume mixing ratio for all species in all model layers of 
less than 10~ 4 to ensure convergence for all atmospheric profiles. Moreover, the height of the cold 
trap is now given by the climate module, and the mixing ratios of all chemical species are adjusted 
according to how much water is evaporated at the surface. 

A main objective of this work is to present simulated spectra for transiting planets, correspond- 
ing to our simulated model scenarios. We therefore use the modelled T-p profiles together with 
the corresponding chemical abundances to calculate spectra for primary and secondary transit ge- 
ometry. The radiative transfer model SQ uIRRL (Schwarzschild Quadrature InfraRed Radiation 
Line-by-line. ISchreier & SchimpflEOOll) was used to calculate synthetic spectra in the infrared 
from 1 .8 to 20 ^m. SQuIRRL was found to be com parable with other radiative transfer codes (e.g. 



von Clarmann et al 



2003 



Melsheimer et al 



20051) . The model takes as input temperature, pres- 



sure, and concentration profiles of the following species from the atmospheric model: KbO, CO2, 
CH 4 , 3 , N 2 0, CO, S0 2 , N0 2 , NO, HC1, HNO3, H 2 2 , H0 2 , and CH 3 C1. H 2 continuum absorp- 
tion corrections are performed and local thermodynamic equilibrium (LTE) is assumed. Absorption 
cross secti ons for all molecules except H2O were computed using the Hitran 2004 line para meter 



data base (IRothman et al]l2005l) . for H2O line data from Hi temp 1995 (IRothman et al 



1995) were 



used. The emission spectra were calculated with a pencil beam looking downwards with a viewing 
angle of 38° to account for the average illumination. Transmission spectra are calculated for dif- 
ferent atmospheric heights (62 layers in total) in limb geometry. The resulting spectra are binned 
for resolutions (R = ^) of R = 5 - 2000, corresponding to instrumental designs for low- and 
medium-resolution spectroscopy. 



2.2. Modelled scenarios 



We study planets with Earth-like chemical composition in the habitable zone of M-dwarf stars 
of types M0 - M7. The spectral appearance of the quiet M-dwarfs is approximated as black 
body emission curves accor ding to effective temperatures and stellar parameters as specified in 
Kaltenegger & Traubl d2009l) (see TableQ}. But these input spectra neglect that many M-dwarf stars 
show strong activity in the UV range, much more than our Sun, for exam ple. The high UVradia- 



tion can have a significant effect on atmospheric chemistry, as discussed in lSegura et al 



(2005). To 



illustrate the effect of high stellar U V flux on the expected spectral appearance of transiting planets, 



4 



H. Rauer et al.: Potential biosignatures in super-Earth atmospheres 
we therefore also u sed a composite spectrum of AD Leo (M4.5V) as central star, as described by 



Segura et al.l (120051) . This high-resolution spectrum is based on observations and a stellar model at- 



mosphere. Observations were taken from satellite dat a (IUE, Inter national Ultraviolet E xplorer) and 
photometry in the visible ( Pettersen & Hawlevl 1989) and near-IR ( Leggett et aljl996 ). Beyond 2.4 



/im, a stellar atmospheric model (NextGen. lHauschildt et al was used to calculate a synthetic 

spectrum. In addition, the Sun was used as a refe rence case f or com parison. The solar spectrum 
was constructed using measurements provided by 



Guevmard 



(2004). The high-resolution spectra 
were binned appropriately for our broadband climate and chemical code. The stellar parameters are 
summarised in Tableland the input spectra used are shown in Fig.Q] Note the significantly higher 
UV flux of AD Leo compared with the quiet M4 and M5 dwarfs shortwards of about 310 nm. 





Fig. 1. Top: Normalised stellar input spectra: Sun (black, solid), AD Leo (dotted), M0 (dashed), 
Ml (dash-dot), M2 (dash-dot-dot-dot), M3 (long dashes), M4 (red), M5 (green), M6 (blue), and 
M7 (magenta). Bottom: Zoom into the photochemically relevant spectral range. 



We checked for the differences in spectral flux between an approximated Planck function and 
real spectral data using the Sun and AD Leo as examples. We found that the differences can be up 
to 20-40 % in the visible and near IR and very large in the far UV. Thus, real stellar input spectra 
should be used whenever possible, especially for detailed chemistry calculations. However, for the 
parametric study performed here the Planck approximation is acceptable for quiet M-dwarfs. 
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Table 1. Stellar parameters used in this study. 



Stellar type 


Teff/K 


R/Ro 


M/M 


Orbital distan 


G2V (Sun) 


5770 


1 


1 


1 


M4.5V (AD Leo) c 


3400 


0.41 


0.45 


0.153 


M7* 


2500 


0.12 


0.09 


0.022 


M6 b 


2600 


0.15 


0.10 


0.030 


M5 b 


2800 


0.20 


0.14 


0.047 


M4* 


3100 


0.26 


0.20 


0.075 


M3 b 


3250 


0.39 


0.36 


0.123 


M2 b 


3400 


0.44 


0.44 


0.152 


Ml b 


3600 


0.49 


0.49 


0.190 


MO 6 


3800 


0.62 


0.60 


0.268 



<a) here the stellar constant equals the solar constant, hen ce the same tot al energy input at TO A is ensured 



Values from 



Kaltenegger & Traubl ( 120091) w Values from lLeggett et all i 19961) 



In our study, planetary gravity was varied from lg to 3g (9.8 to 29.4 m s ). Hence, our l g 



and 3g cases correspond to 1 and 1 1 .4 Earth masses and 1 and 1 .95 Earth radii (ISotir 



respectively. The surface pressure was kept fixed at 1 bar. As in previous studies (e.g. 



et al.H2007) 



2003 



Segura et al 



Segura et al.1 120051) . we set the model parameters surface albedo and fluxes of chemical 



species in a way that the mean surface temperature of 288 K and observed surface concentra- 
tions of the chemical species H2, CO, CH4, N2O and CH3CI are reproduced for present Earth. For 
the other model scenarios we kept the surface albedo and fluxes constant and calculated the atmo- 
spheric concentrations and surface temperature. However, in contrast to previous studies, we then 
placed the modelled extra-solar planets at an orbital distance from their central star where their 
total stellar energy input equals that of 1 solar constant (1366 W/m 2 ). We chose this approach so 
that habitable surface temperatures were not enforced and allowed the temperature to adjust in a 
consistent way to the different spectral energy distribution of the M-dwarf stars. Thus, we separated 
in our investigation the effects of the different spectral energy distributions from variations of the 
total incoming flux. We will show in Sect. 4 that the resulting surface temperatures are still in the 
habitable range. 

At present, the atmosphere of a terrestrial extrasolar planet has not been detected directly. 
Assumptions on their composition are therefore more or less arbitrary, and are based on our knowl- 
edge of the Solar System. In terms of searching for biomarkers, Earth is the only reference case 
we know of a planet with a biosphere. In this study, we therefore assume the same initial modern 
Earth atmospheric composition and biomass emissions for the simulated extra-solar planets. The 
final atmospheric composition for the star-planet-scenarios studied was then consistently calcu- 
lated according to the coupled climate-chemistry model. Although the assumption of an Earth-like 
biosphere for planets orbiting stars that are quite different from our Sun is probably unrealistic, 
it nevertheless serves as a first step towards understanding the effects of different environmental 
factors on habitable planets and their atmospheres. 

We used specifically the following surface boundary conditions to reproduce mean Earth con- 
ditions in our Sun-Earth reference run. Molecular hydrogen was removed at the surface with a 
constant deposition velocity of 7.7xl0~ 4 cm s _1 . Biomarkers (N2O and CH3CI) and related com- 
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pounds (CO an d CH4) were emitted w ith constant upwards fluxes per surface area in a procedure 
outlined e.g. in 



Grenfell et al 



(12007bl) . For example a CH4 flux of 531 Tg/year is needed to repro- 
duce modern Earth CH4 surface volume mixing ratios and is kept constant for all other simulated 
scenarios. This approach is different to the one employed bv lSegura et al.l (120051) . who kept the CH4 
surface mixing ratio for the quiet M-star runs constant. Furthermore we did not need to change the 
C O surface boundary condition from a constant flux to a constant deposition velocity as was done 
Segura et al ( 2005 ) for their qu iet M-star scenarios. As a result, our surface fluxes differ largely 



from those of 



Segura et al.l (12005b but compare better with the observed values on Earth. Again 



we chose our approach to vary input parameters step-by-step, not to mix the influence of different 
effects in our analysis. 

To validate our approach, we performed S NR calculations for several m odern Earth transmis- 
sion features, which were also performed by iKaltenegger & Traubl (120091) (their Table 2, single 
transit events). When we strictly followed their approach (fixed the atmospheric T-p profiles and 
composition to modern Earth), we were able to reproduce their values to within 2%. A direct com- 
parison of the model results with observations of the Earth transmission spectra by e.g. lPallel (120101) 



is unfortunately not possible. The measured spectrum is shortwards of 2.4 jum and thus overlaps 
only little with our modelling range. Furthermore, because our model is stationary and cloud-free 
a comparison wi th real data of a clo udy Earth is difficult and no direct agreement is expected. A 
comparison with lSegura et al.l (120051) is discussed in more detail in Sect. 5. 



3. SNR calculations 

We here consider emission and transmission spectroscopy in the near- to mid-IR wavelength range. 
We account for photon noise only to give an upper limit on the detection feasibility. Additional 
noise sources, e.g. instrumental or astrophysical, like zodiacal light, will further decrease the signal- 
to-noise ratio (SNR) computed here. This section describes the method used to calculate the respec- 
tive SNR. 



3.1. Emission 

The planetary emission spectrum, F p , is obtained by combining fluxes (with equal integration 
times) observed during (F\) and outside (F2) secondary eclipse 

F P =F 2 -F U (1) 

where F] = F s and Fo = F s + F p , and F s represents the stellar flux. The SNR for the planetary 
spectrum is therefore 

SNR p = , Fp (2) 



+ 07, 



Here, cr s and cr p denote the stellar and planetary flux standard deviation, respectively. Assuming a 
photon-limited instrument (hence, only Poisson noise from star and planet) and F s » F p , we then 
reduce (2) to 

Fp F s C E 

SNR p = — J— = C E ■ = SNR, ■ (3) 
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with SNR, the stellar SNR and Ce the contrast between planetary and stellar signal in the secondary 
eclipse emission spectrum. The resulting SNR is conservative compared with the continuum flux 
region. 

Ce is calculated from the integrated fluxes of a hemisphere as 

F„ R\ I P 

where Rp and R$ are planetary and stellar radius, respectively, and Ip and Is are the spectral en- 
ergy fluxes (in W m 2 //rrT 1 ) from planet and star, respectively. Here, the st ar is taken as a black- 



body with a fixed effective temperature, consistent with previous studies (e.g. jKaltenegger & Traub 
20091) . The observable stellar signal S s (in J junT 1 ) is calculated by 



R 

S s = -p'Is-A't = S« U 'A-t, (5) 

where d is the distance to the star, A is the telescope area (assuming a 6.5 m JWST-like aperture 
in this case) and t the integration time. To convert into the number of photons, we divide S s by 
N - hj (h Planck constant, c speed of light, A wavelength). From there, we can calculate the 
spectrally integrated stellar signal S j (number of photons) 

7 ~ AT ' ¥' (6) 
with Ac the central wavelength and R the specified spectral resolution. 
S det, the signal received from the star by the detector, is calculated as 

S de , = 5/ • q, (7) 
where q is a factor to acc ount for a non-ideal detecto r and the overall instrument throughput (as- 



sumed for JWST: g =0.15 lKaltenegger & Traubll2009l) 



The signal-to-noise ratio of the star, SNR S , is calculated without taking into account instrumen- 
tal or read-out noise. 



SNR S = ^= V^= Js* t -A-t--2£--q. (8) 
cr s V N ■ R 

The SNR values are calculated for integration times equal to the transit times of considered 

scenarios (f = To)- The duration, To, of a transit is calculated by (assuming that Rp <K R$, i = 90°) 



. (Rs 
To - — -2 arcsin 
2tt 



where P is the orbital period and a the planetary orbital distance. 
From Eq.|3]we then calculate the planetary SNR (SNR P ). 

Integration times f, for a certain specified SNR,,. are obtained from 



/ SNR, . 

^MsnrT 1 (10) 



using the yi scaling of SNR with integration time. 
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Because the SNR scales with VF and therefore d~\ the integration times calculated for a 
specific distance d\ can be easily scaled to any distance d?. 

t{d 2 ) = t{d{)-{^\ . (11) 

3.2. Transmission 

The transmission contrast is, to a first order approximation, the photometric transit depth, i.e. 

Ct = -f , (12) 
R s 

which is much higher than Ce from Eq. [4] 

In our transmission calculations, we simulated 62 adjacent tangent height rays through the 
atmosphere, with a diameter corresponding to the altitude layer difference in the chemical module. 
The central altitudes of the beams correspond to the layer altitudes of our model atmosphere. From 
there, the effective height h(A) of the atmosphere is calculated as 

A(A)=2(l-r<(A))M,-, (13) 

i 

where T, is the transmission for ray i at wavelength A and Ah, the altitude difference between 
consecutive tangent heights. The additional transit depth, fa, provided by the atmosphere is there- 
fore obtained by 

^±*C^_l). (14) 



The overall transmission, T, of a beam at a given wavelength A is calculated by 



r^^^i, (15) 



where H is the total height of the considered atmosphere. Because our beams cross the atmosphere 
in equidistant grid points, T is the simple arithmetic mean of all T,. 

The SNR of a spectral feature in transmission is then calculated by (analogous to Eq. [3} 

SNRj- = SNR S ■ -4=. (16) 
V2 



Note the additional factor -4= in our SNR compared with iKaltenegger & Traubl (12009b . The 
factor results from our approach that considers the planet signal as resulting from the difference of 
the measured stellar flux in and out of transit. 

Required integration times or integration times vs. distance are then obtained by the same scal- 
ing as in Eqs.fTOlandfTTI 



4. Results 

4.1. Atmospheric profiles 
4.1.1. Temperature profiles 

Figure [2] shows the temperature profiles for two different planet gravity values (lg, 3g) with solar 
radiation in both cases. The temperature profile for the modern Earth control run shows a convec- 
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180 200 220 240 260 280 300 
T/K 

Fig. 2. Influence of gravity: Earth-sized (solid) vs. 3g super-Earth planet (dashed) around the Sun. 

tive troposphere, and the temperature inversion and stratospheric temperature maximum caused by 
ozone heating. 

The surface temperature decreases for the 3g scenario. This is because we held the surface 
pressure constant, which consequently lowered the vertical column mass. Less mass means less 
absorption and therefore less greenhouse effect, hence lower surface temperatures. Furthermore, 
effective ozone heating occurs at higher pressures, p, and the temperature maximum in the strato- 
sphere accordingly decreases to higher p when increasing gravity. Overall increased stratospheric 
ozone levels (see Sect. 4.1.2) lead to enhanced heating in these parts of the atmosphere. 




Fig. 3. Influence of stellar spectrum: Earth-sized planets around the Sun (black, solid), AD Leo 
(dotted), M0 (dashed), Ml (dash-dot), M2 (dash-dot-dot-dot), M3 (long dashes), M4 (red), M5 
(green), M6 (blue), and M7-type stars (magenta). 

Figure[3]shows the influence of the stellar type on the temperature profile for the Earth control 
case as well as lg planets around AD Leo and the M0 to M7 stars from Table Q] Because unlike 
in previous approaches, we did not fix the surface temperature of our modelled planets by different 
stellar insolations at TOA, it is interesting to note that although Earth-like planets around M-dwarfs 
show somewhat increased surface temperatures, they always remain well within the habitable range 



10 



H. Rauer et al.: Potential biosignatures in super-Earth atmospheres 



for a total received energy flux of 1 solar constant. However, surface temperatures for the lg sce- 
narios are up to 20 K higher than for the Earth-Sun case. For the super-Earth scenarios similar 
changes occur. Here, surface temperatures are higher for the M-type central stars than for the Sun, 
and increase from MO to M4-type stars (not shown). The increase occurs because Rayleigh scatter- 
ing contributes less to the planetary albedo owing to the /Independence of the Rayleigh scattering 
coefficient. In addition, we find enhanced greenhouse warming (see below). Overall, however, the 
differences in spectral energy distribution alone have only a moderate effect on mean surface tem- 
peratures in the context of habitability in general. They may nevertheless significantly alter the 
planetary environment when studying the detailed atmospheric condi tions (for exampl e in a 3D 
approach that includes dynamical effects such as tidal-locking (see e.g. Joshi et alll997 )). 



The main di fferences between th e Sun and the M-star planets are found in the stratosphere, as 
already noted bv lSegura et al.l (120051) . The atmospheres of the modelled M-star planets do not show 
a stratospheric temperature inversion up to type M5 (M4) for the Earth-sized (super-Earths) planets. 
This is caused by the different stellar energy flux distribution of M-dwarfs compared with the Sun. 
The reduced UV flux longwards of about 250 nm results in less heating by ozo ne and therefore a 



reduced stratospheric inversion. We confirm this finding by ISegura et al 



(2005) for the active star 



AD Leo. In addition, our study also includes model runs for several quiet M-dwarfs with otherwise 
the same parameters for their orbiting planets. As can be seen in Fig.[3j the planetary stratospheric 
temperatures increase from hot (M0) to cooler M central stars (M7), even building up a strato- 
spheric inversion again for M6 and M7 stars. This is caused by increasing amounts of methane 
and water in the mid-atmosphere (see following section) which lead to enhanced absorption and 
therefore heating of the stratosphere. 



As already noted by ISegura et al 



von Paris et al 



(2008), the RRTM radiation 



(120031) and 

scheme as applied in the model has been constructed for atmospheres that do not deviate too muc h 
from Earth conditions both in terms of temperature structure and composition ( Mlawer et al.l 1997 ). 
When applied to atmospheres that significantly deviate from modern Earth, it is important to check 
that the atmospheric T-p-profiles do not leave the validity range of RRTM. If this is the case, the 
resulting effect on the heating and cooling rates must be estimated. 

We compared the computed T-p-profiles for the atmosphere considered here to the validity 
range of RRTM. We found that T-p-profiles for most M- dwarf central stars lea ve the validity range 
significantly only in the upper atmosphere. As shown bv lvon Paris et al.l (120081) . variations in upper 



stratospheric optical depths have a negligible influence on lower atmospheric temperature profiles. 
Most of the emission spectra originate in pressures higher than 0. 1 bar, except for planets around 
very cool M dwarfs. For a given chemical composition, the transmission spectra are nearly inde- 
pendent of T-p structure (in contrast to the emission spectra), hence uncertainties in the IR radiative 
transfer of the climate model are unlikely to significantly influence the transmission spectra. To in- 
vestigate whether these temperature errors might in turn influence the chemical composition of the 
atmosphere by means of e.g. temperature-dependent reaction rates, additional tests were performed 
in which the stratospheric temperatures were arbitrarily changed by +10K. Observed changes in 
concentrations and column amounts did not exceed 5-20% depending on species, with methane 
showing the largest effect. However, no difference could be discerned in the spectra. Hence, the 
influence of errors in stratospheric optical depths on the planetary spectra is likely to be small. 
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A comparison of the thermal fluxes of RRTM with thermal fluxes computed by the line-by-line 
radiative transfer code SQuIRRL showed an agreement of the total outgoing fluxes within 5.5%. 

The surface temperature decreases again for the lg planet that orbits the M7 star. In view 
of the above discussion on model validation, we consider the M7 case to be at the limit of the 
application range. We did not study cooler central stars here, because the results would be out of 
the application range over most of the atmosphere. For the same reason, we also exclude the M6 
and M7 stars from our further analysis of 3g super-Earth scenarios. Nevertheless, the M7 case 
indicates an interesting trend. Increased absorption of stellar radiation in mid-atmosphere levels 
owing to massively increased H2O and CH4 concentrations results in less radiation reaching the 
ground. Whether this trend would continue for planets that orbit even cooler M dwarf stars would 
need to be further investigated with a model adapted to these conditions. 



4.1.2. Chemical profiles 



Ozone, nitrous oxide, methyl chloride, methane, and water are the major biomarker molecules and 
related compounds. In this section we give a broad overview of the concentrations of O3, N2O, 
CH3CI, CH4 and H2O. A detailed analysis of th e chemical processes for the model runs presented 
here will be given in a forthcoming publication (IGrenfell et alJ (I2010I) . hereafter paper II). 
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Fig. 4. Chemical profiles for an Earth-sized planet (solid) and a 

3g super-Earth (dotted) around the Sun. 



Chemical profiles for 1g and 3g scenarios 

Because of the decrease in total column mass for the 3g planetary scenario compared with the 
lg case (because of the fixed surface pressure of 1 bar), the decrease in greenhouse effect (GHE) 
resulted in lower tropospheric temperatures. This in turn led to a small lowering in the tropospheric 
H2O number densities, as can be seen in the upper left panel of Figure|4] Although the impact on the 
number density is small, the impact on the H2O column is large (it is five times lower for the 3g case 
than for the lg scenario) owing to the overall decrease in column mass. The decrease in the H2O 
column led to a lowering in OH to a surface number density of 1.24T0 6 molec. cirT 3 compared 
with 1 .3 1 TO 6 molec. cirT 3 for the lg case. Because OH is a major sink for CH4, its number density 
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and column increased. While the surface CH4 number density increased from 3.9T0 13 molec. cm -3 
(lg) to 1.7-10 molec. cm -3 (3g) (see upper right panel in Figure|4), the column increased by 44% 
only, again owing to the lower total column mass for the 3g scenario. 

The O3 number density increased for the 3g planetary scenario (see lower left panel in Figure 
4) in the lower to mid-atmosphere compared with the lg case. The overall smaller column amount 
at 3g led to a small strengthening of the Chapman source at these levels (for more details see paper 
II). The main ozone photochemical features typical of the Earth were maintained at lg and 3g, e.g. 
in both cases HOx and NOx catalytic cycles were important in the upper and lower atmospheres 
respectively, whereas the smog cycle was important for ozone production in the troposphere for 
both lg and 3g. Compared with changing the stellar type (see below), the effects of increasing 
gravity are rather small, however. Nevertheless, a more detailed study of the chemical effects in 
simulated super-Earth planet atmospheres is given in paper II. 
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Fig. 5. Influence of central star type on chemical atmospheric profiles of an Earth-sized planet: Sun 
(solid), MO (long dashed), AD Leo (dotted), M5 (dot-dashed), M7 (triple-dot-dashed). 



Chemical profiles for different stellar types 

The response of biomarkers and important greenhouse gases to different central stars is shown 
in Fig. [5] for some examples. Increasing the M-dwarf stellar type and consequently decreasing 
the amount of UV radiation (except for the active star AD Leo) has a significant effect on the 
photochemical responses in the atmosphere. We briefly discuss the main photochemical responses. 
For more details we refer to the analysis of the chemical processes made in paper II. 

We first note a strong increase in planetary stratospheric water content when moving from the 
Sun via MO stars to our coolest central star, M7. Furthermore, methane significantly increases for 
cooler central stars. One cause for increased stratospheric water is the increase in tropospheric wa- 
ter abundance owing to higher surface temperatures. Another reason is the reduced UV flux of quiet 
M dwarfs compared with the Sun, which leads to less photolytic destruction (see Fig. [TJ. As a re- 
sult, the abundance of OH is also reduced, which can lead to reduced CH4 destruction. On the other 
hand, increased methane abundances can lead to faster CH4 oxidation and therefore to increased 
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Fig. 6. Contrast spectra for an Earth-sized planet (1 g) around the Sun (black), AD Leo (red), MO 
(green), M5 (blue), and M7 (magenta). 



stratospheric water production, as suggested bv lSegura et alj d2005l) . Detailed investigations show 
that other catalytic cycles also play a role. This analysis is presented in paper II. 



Segura et al.l (12005b pointed out that the photolytic destruction of methane in Earth-like, high- 



oxygen planets is slow, and in combination with drastically reduced sinks by OH this could lead to a 
"methane runaw ay". Possible m ethane runaway feedbacks are quite well investigated for the Earth's 
atmosphere (e.g. 



Prather 



19961) . The main feedback is positive, where more methane leads to less 



OH, therefore slower methane in-situ oxidation, which favours methane even more. There is also an 
opposing, negative feedback on Earth, e.g. methane warming leads to a damper atmosphere, which 
favours an increase in OH. There exist potentially complex interactions e.g. with NOx emissions, 
where increased NOx generally favours fast parti tioning of HO2 i nto O H (via: NO + HO2 — > OH 
+ NO2) which in turn affects the CH4 runaway. 



Lelieveld et al 



(2002) provid ed an overview of 



these effects. In model simulations as presented here and in lSegura et al.l ( 12005b . methane runaway 
can occur as an unphysical scenario for quiet M-dwarfs if the assumed methane fluxes exceed 
the chemical sinks and loss processes in the model. Our methane fluxes are lower than those of 
Segura et al.1 (120051) . and this methane runaway does not occur for our scenarios. 

Fig. [5] shows the ozone profiles for the Earth-sized planet scenarios. The ozone profile of the 
Earth-Sun reference run and the planets that orbit AD Leo are similar, whereas ozone decreases 
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for cooler M-dwarfs and is significantly reduced in the M4 to M7 cases. Ozone abundances are 
controlled by the balance between photochemical sources via the established Chapman and smog 
mechanisms and the catalytic O3 sinks, which depend on concentrations of e.g. hydrogen and 
nitrogen-oxides. A key result of the analysis of chemical processes in paper II is that ozone pro- 
duction changes from being dominated by the Chapman to the smog mechanism with increasing 
stellar type, thus with decreasing effective temperature and UV flux. The weaker UV flux from the 
cooler M dwarfs leads to slower photolysis of the O2 molecule, which is required by the Chapman 
mechanism. Smog production is still slow in comparison. This explains the decline of ozone for 
cool, quiet M-dwarf central stars. The simulated planet that orbit the active star AD Leo though 
shows a pronounced ozone layer similar to the Earth owing to its higher UV flux compared with 
the quiet M4 and M 5 stars. The prominent ozone layer for active M dwarfs was also found by 



Segura et al.l (12005b - More details on the changing ozone production mechanisms for cool central 



stars are explained in paper II. 

The changing stellar flux of M dwarfs with increasing M-star class also leads to a strong in- 
crease in N2O densities (Fig. [5]). This is again caused by reduced photolysis rates due to the lower 
UV flux. Again, the effect is smaller for active M-dwarfs like AD Leo, where the photolysis rates 
are much higher. 

W e also investigated t he abundance of the biomarker molecule CH3CI (not shown here), be- 
cause ISegura et al.l (120051) showed this to potentially cause significant ab sorption bands for high 



concentrations and active M-dwarf stars. In contrast to lSegura et al.l ( 12005b . however, our required 
surface CH3CI flux to reproduce mean Earth conditions is a factor of 2 lower, which also better 
agrees with the observed estimates. We find an increase in CH3CI abundances for cooler M dwarfs 
owing to the reduced UV flux, similar to the increase in methane. However, in all cases studied, no 
significant spectral absorption could be seen in the emission spectra, whereas small absorption at 
13.7 /mi for the coolest central star (M7) can be found in the transmission spectrum (see Fig. [Tot. 
This molecule was therefore not investigated in more detail in this work. 



4.2. Emission spectroscopy 
4.2.1. Emission spectra 

The IR emission spectra of extra-solar planets can be measured for transiting planets by observing 
the secondary eclipse. The planetary signal is derived from the difference of flux during and out- 
of-eclipse. Below, we show the expected flux contrast (planet/stellar flux) for computed planetary 
IR emission spectra. 

The modern Earth reference contrast spectrum is shown in Fig . [6] (black line). The spectrum 
follows a blackbody spectrum with absorption bands of, e.g., water (2.7 /mi, the broad 6.3 //m band 
and the rotation band longwards of 15 /mi), methane (at 3.3 and 7.7 //m), ozone (9.6 /mi) or carbon 
dioxide (2.7, 4.3 and 15 /mi). The contrast between stellar and planetary flux is very low, reaching 
10" 6 in the mid-IR. Because of the rapid decrease of the planetary Planck function towards shorter 
wavelengths, the contrast falls below 10~ 10 in the near-IR. 

Figure [6] also compares the spectra of an Earth-like planet (lg) when changing the central star 
(Sun, AD Leo, M0, M5, M7). The overall contrast is improved for terrestrial planets around M 
dwarfs, as expected. Because of the much lower effective temperatures and the smaller radius of 
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Fig. 7. Contrast spectrum of the 9.6 fim ozone band and the 15 yum carbon dioxide band for an Earth- 
sized planet around the Sun (black), AD Leo (red), MO (green), M5 (blue), and M7 (magenta). Blue 
diamonds: binned to R=20, red: smoothed spectra at the same resolution. 



M-dwarf stars, the contrast increases by about two orders of magnitude and eventually reaches 
about 10 in the mid-IR, which illustrates the improved potential of detecting planets around M- 
type dwarf stars. 

Concerning molecular absorption bands, noticeable differences to the Earth-Sun system are 
found in the 2.3, 3.3 and 7.7 fim methane bands and the 3.8 and 4.5 jum band of nitrous oxide. 
These molecules show higher concentrations in the atmosphere of M dwarf planets compared with 
Earth (see Fig. [5J, therefore we would expect their spectral signatures to be more pronounced. 
Additionally, the 4.7 pm band of CO becomes more prominent owing to its higher concentrations 
in the atmosphere (not shown). This can indeed be observed for the planets that orbit the quiet 
M0-M3 dwarfs stars. However, the absorption bands of N2O, CH4 become weaker for cooler M 
stars and almost completely disappear for an Earth-sized planet around the M7 star, despite the 
greatly enhanced abundances. This is caused by the temperature structure of these planets (see 
Fig. [3}. The emission in the CH4 and N2O bands originates in the middle atmosphere because 
of the large concentrations of these gases, i.e. the atmosphere becomes transparent only at much 
higher altitudes than for Earth around the Sun. At these altitudes, the temperature is almost as 
high as the surface temperature, hence the contrast between continuum (transparent atmosphere) 
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Fig. 8. Emission spectra: Integration time in hours needed for SNR=3 as a function of stellar dis- 
tance of the 9.6 fim (dotted) ozone and the 15 fim (solid) carbon dioxide spectral feature at R — 20 
for an Earth (black) and a 3g super-Earth (red) around AD Leo (top) and an M5 dwarf central star 
(bottom). Telescope configuration: JWST. 



and absorption band is reduced. The absorption cannot be discerned from the spectrum any longer. 
Thus, the atmospheric composition may be difficult to detect via emission spectroscopy in planets 
orbiting quiet and very cool M dwarfs, despite a more favourable planet/star flux contrast ratio. 

Another remarkable difference occurs for the 4.3 and 15 fim fundamental bands of carbon 
dioxide. The absence of a temperature inversion and the much cold er stratosphere for th e MO to 
M5-type stars (see Fig.0 have a noticeable effect, as already noted bv lSegura et al.l (120051) . The 4.3 



fim band absorption is stronger for the MO and AD Leo cases, and in the 15 fim band the emission 
peak is absent. However, as already discussed, our model runs for quiet M stars show that for cool 
stars, e.g. the M7 case, the emission spectrum originates mainly in the upper atmosphere for almost 
the entire spectrum (not only the CO2 fundamentals), and therefore molecular absorption features 
are weak. 

Figure [7] shows the 15 fim carbon dioxide and the 9.6 fim ozone band at high and low wave- 
length resolutions. The CO2 band is clearly visible already at a resolution of R = 20, although 
the emission peak close to the line centre is not resolved. This emission peak is an indicator of 
the atmospheric temperature inversion, hence a potentially important spectral feature. To resolve 
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Fig. 9. Effective tangent heights: Earth and 3g super-Earth (red) around Sun. 



the emission peak, however, at least R « 200 is necessary. The same conclusion holds for other 
interesting very narrow features. The ozone absorption at 9.6 fim is still visible at R — 20. For a 
detection of the molecule, low resolution seems to be sufficient. We also performed simulations 
for even lower resolution {R = 5). In this extreme, however, ozone is barely visible against the 
surrounding continuum and becomes difficult to detect. 



4.2.2. Signal-to-noise-ratios 

We calculated the expected SNR for different planet-star c onfigurations. The chosen tel escope aper 



ture is as for JWST with an assumed efficiency of 0.15 (Kaltenegger & Traub 2009). Integration 
times correspond to a single transit event (i.e., 13 hours for super-Earth around Sun and 3 hours for 
super-Earth around AD Leo). For the Sun-Earth case, SNR values (not shown) are always below 
0.15. Earth-like planets around M stars also show low SNR values. Table[2]summarises the achiev- 
able SNR values for a 3g super-Earth around the Sun, AD Leo, M0, and M5. As can be seen from 
Table|2j the single transit SNR is below unity for all absorptions, except for low resolution and late 
M dwarf types. 

Table [2] shows that emission spectroscopy of terrestrial planets in the HZ provides very low 
SNR for a single transit duration exposure for the 3.3 fim methane, the 4.3 fim carbon dioxide, the 
6.3 fim water and the 7.7 fim methane band. However, for the 9.6 fim band of ozone and the 15 /an 
band of carbon dioxide, SNR values are more encouraging at low resolution. Nevertheless, several 
transits need to be added to reach realistic SNR. 

In Fig. [8] we show the integration times needed to obtain a SNR of 3 for R — 20 as a function 
of stellar distance. The planetary scenarios are the lg Earth-like and 3g super-Earth planet around 
AD Leo and an M5 dwarf. As examples, we chose the ozone band and the broad carbon dioxide 
band at 15 fim. The integration time is indicated in units of transit time on the right y-axis. The 
figure indicates that for 3g super-Earths a SNR of 3 in the 15 fim CO2 band could be achieved in a 
single transit for the cool M5 dwarf closer than about 10 pc. Higher resolution is probably required 
for more detailed characterisation, e.g. for the detection of a temperature inversion. 
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Table 2. SNR values for emission spectra of a 3g super-Earth at 10 pc" 
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Integration time is a single transit duration. Telescope configuration: JWST. 



4.3. Transmission spectroscopy 
4.3.1. Transmission spectra 

During primary eclipse the light of the central star passes through the atmosphere of the transit- 
ing planet and absorption by atmospheric molecules will change the apparent radius of the planet 
according to their scale height. Generally, transits are observed in the optical wavelength range 
because of the higher stellar flux and therefore better SNR. However, M-type stellar radiation 
is shifted somewhat to IR wavelengths and, furthermore, many interesting biomarker or related 
molecules show absorption in the IR. Below we therefore investigate the signal for primary transit 
spectroscopy in the near- to mid-IR range. 

The wavelength-dependent effective tangent height of the planetary atmosphere blocks the stel- 
lar light, hence increases the apparent radius of the planet. The effective height /; (see Eq. [T3ll is 
shown in Fig. [9] for an Earth-sized planet (black line) and a super-Earth (red) around the Sun. h 
is mainly below 20 km, except for the strong methane 3.3 fim band, the ozone 9.6 /urn band and 
the three 2.7, 4.3 and 15 fim bands of carbon dioxide in the case of the Earth. The super-Earth 
scenarios, however, may not be best suited for transmission spectroscopy, because the effective 
tangent height is of the order of the atmospheric scale height, which varies as g - 1 . From E q.fT4l the 
planetary signal decreases as g 0A when using the mass-radius-relationship from 
This effect is illustrated in Fig. [9] 



Sotin et al 



(2007). 
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Fig. 10. Comparison of atmospheric transmission of a lg Earth-like planet around the Sun (black), 
AD Leo (red), MO (green), M5 (blue), and M7 (magenta) dwarf stars. 



The effect of changing the central star for the lg scenarios is shown in Fig.[l0]for various M- 
dwarf type central stars. The relative transmission spectrum (from Eq.Qj]) of the Earth is shown in 
black. Absorption from water, ozone, carbon dioxide, and methane are clearly visible. Furthermore, 
the nitrous oxide 4.5 fim band and 17 fim band as well as a HNO3 band at 1 1.5 /mi can be distin- 
guished. Some interesting differences in the spectra are seen for M dwarfs. The 15 fim band of 
CO2 changes only by a small amount from the Sun up to the M7 central star. However, the ni- 
trous oxide, methane and water bands continuously increase from earlier (hotter) to later (cooler) 
M dwarfs. We note that N2O absorption can be as strong as the main CO2 bands for the simulated 
planet around the M7 star. In this case, the nitrous oxide would be a strong biomarker signal. The 
HNO3 absorption basically disappears for AD Leo and is weak for the central MO star case, but 
strongly increases up to the quiet M7 star. 

Ozone is strongest in AD Leo and increases somewhat from MO to M3 for the quiet central 
stars. The ozone signal then, however, decreases again for planets around cooler central stars due 
to the reduced overall O3 column. Thus, although atmospheric absorption signals in transmission 
may be easier to detect for planets orbiting cooler stars for most species, this may not be the case 
for ozone. 
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Table 3. SNR for transmission spectra of an Earth-like planet at 10 pc" 
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1.610 1 


7.8-10 2 




C0 2 


4.3 


1.24 


9.M0-' 


2.0-10-' 


9.9- 10" 2 




CH 4 


7.7 


8.5-10-' 


4.9-10-' 


1.1-10"' 


5.3-10- 2 




3 


9.6 


5.5-10-' 


4.3-10-' 


8.510 2 


3.3-10- 2 




C0 2 


15.0 


5.4- 10-' 


3.110 1 


6.810 2 


3.3-10- 2 


MO 


CH4 


3.3 


8.3-10-' 


5.4-10-' 


1.1-10"' 


4.9- 10- 2 




C0 2 


4.3 


9.6- 10"' 


8.110 1 


1.8-10" 1 


8.4- 10- 2 




CH 4 


7.7 


5.2-10-' 


3.3-10' 


7.110 2 


3.5-10- 2 




3 


9.6 


4.2-10-' 


3.6-10-' 


7.110 2 


3.7- 10- 2 




C0 2 


15.0 


4.5-10-' 


2.8-10-' 


5.9-10- 2 


2.9- 10- 2 


M5 


ca, 


3.3 


1.94 


1.19 


2.1-10" 1 


1.0-10- 1 




co 2 


4.3 


1.86 


1.20 


2.6- 10"' 


1.3-10" 1 




CH 4 


7.7 


1.27 


7.4- 10-' 


1.610"' 


8.0- 10- 2 




o 3 


9.6 


6.2- 10-' 


4.4- 10-' 


8.4- 10- 2 


4.4- 10- 2 




C0 2 


15.0 


7.2- 10"' 


4.2- 10-' 


9.2- 10- 2 


4.5-10- 2 


M7 


CH 4 


3.3 


2.66 


1.59 


2.9-10' 


1.4-10-' 




C0 2 


4.3 


2.47 


1.40 


3.0-10-' 


1.5-10-' 




CiU 


7.7 


1.78 


1.02 


2.2-10-' 


1.1-10-' 




3 


9.6 


7.7- 10-' 


4.7-10-' 


9.110 2 


4.7- 10- 2 




C0 2 


15.0 


9.0- 10"' 


5.0-10"' 


1.1-10 ' 


5.3-10- 2 



Integration time is a single transit duration. Telescope configuration: JWST. 



4.3.2. Signal-to-noise-ratios 

Table [3] summarises the achievable SNR for different absorption bands for transmission spec- 
troscopy during a single transit for Earth-sized planets around the Sun and M dwarfs. In contrast 
to emission spectroscopy, transmission spectroscopy does not probe the temperature structure of 
an atmosphere, but the chemical composition. Accordingly high resolution is not explicitly neces- 
sary. Unfortunately, at mid-IR wavelengths (e.g. 15 yum), the stellar SNR is relatively low, which 
makes detection harder than at shorter wavelengths (e.g., ozone at 9.6 yum or methane at 3.3 /mi). 
For near-IR absorption bands of CH4 and CO2 a SNR>1 can be achieved in a single transit, thus 
providing the possibility to achieve reasonable SNR if several transits can be co-added. Figure [TT] 
illustrates this concept. It shows integration times as a function of required distance to achieve a 
SNR of 3. 

Figure QT| implies that for close stars (closer than 5 pc) a JWST aperture-like telescope would 
be capable of performing low-resolution transmission spectroscopy for Earth-like planets on all 
bands for the simulated M5 dwarfs. Earlier types of M dwarfs, however, provide more challenging 
detection conditions. To cross-check our results and identify the differences of the approach used 
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Fig. 11. Transmission: Integration time needed for SNR=3 as a function of stellar distance of the 
3.3 fim methane (dash-dot-dot-dot), 4.3 /mi carbon dioxide (dash-dot), 7.7 fim methane (dashed), 
9.6 fim (dotted) ozone and the 15 fim (solid) carbon dioxide spectral feature at R = 20 for an Earth 
(black) and a 3g super-Earth (red) around AD Leo (top) and M5 (bottom). Telescope configuration 
is JWST 



here, we compared our approach to that used by iKaltenegger & Traubl (120091) . The difference is 



that we calculated the atmospheric co ncentration profiles consistent with different M dwarf spectra, 
whereas lKaltenegger & Traubl (120091) used a fixed modern-Earth reference profile. Both approaches 



yield quite similar results for the CO2 15 fim band. However, for e.g. the methane and ozone 
absorption we find significant differences. This shows that a consistent approach is necessary to 
take into account the effect of stellar spectral energy distribution on atmospheric T-p profiles and 
biomarker mixing ratios. 



5. Summary and conclusions 

We investigated the effect of the stellar energy flux distribution of M-dwarf stars on atmospheric 
profiles of Earth-like and super-Earth planets. In addition, IR emission and transmission spectra 
and SNRs were computed to investigate how this affects the spectral appearance and detectability 
of absorption bands. We studied planets with lg and 3g (1-2 Earth radii) because this approximately 
corresponds to the range from Earth to super-Earth planets. The planets were placed at a distance 
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from their host star where the stellar energy input equals the energy of the Earth received by the 
Sun. To illustrate the detection capabilities for planets with biospheres, we arbitrarily assumed 
planets with modern Earth atmospheric composition as an initial starting value. However, the actual 
atmospheric composition was then re-calculated according to the incident stellar spectral flux to 
obtain a consistent model of the climatic conditions and composition of the atmospheres for these 
planets. 

Our results suggest that the different spectral energy distribution of M dwarfs has a potentially 
large effect on the atmospheric T-p-profiles and on the chemical profiles of H2O, CH4, O3, N2O, 
and HNO3. In particular the temperature inversion i n the mid-atmosphe re disappears for planets 
around M dwarfs (Fig. |3}- This was already found bv lSegura et al.l (12005b for planets around active 
M-dwarf stars. For quiet stars with low UV flux we find, however, that an inversion can build 
up again in the planet's atmosphere owing to enhanced concentrations of water and methane in the 
mid-atmosphere. In contrast to previous approaches, we did not fix the surface temperature to 288K 
through stellar flux adjustments in our simulations. Surface temperatures can be up to 20 K higher 
for planets around M dwarfs for the same total stellar energy input, but the simulated planets always 
retain their habitable conditions. The effect on atmospheric structure and composition appears to 
be weaker when we varied gravity from l-3g compared with varying the stellar input spectra (Fig. 

m. 

We found systematic increases in the atmospheric abundances for water, methane, and nitrous 
oxide owing to the reduced photolysis rates associated with M-dwarf stars compared with the Sun. 
Abundances generally increase from hot to cool M dwarf stars. For active M dwarfs, like AD Leo, 
photolysis heating rates are higher, and thus their related chemical processes are stronger than 
for the quiet M dwarfs of a similar type. Interestingly, planetary ozone abundances decrease for 
very cool, quiet cen tral stars. This is cau sed by low photolysis rates which reduce the Chapman 



ozone production. In lGrenfell et al.l (120101) (paper II) we present a detailed analysis of the chemical 



responses of these atmospheres. 

The main goal of this work is to present simulated spectra and SNR estimates for the modelled 
planetary systems. We therefore discussed simulated spectra for secondary and primary transit ge- 
ometries. During secondary eclipse, the measured emission fluxes depend not only on the composi- 
tion of the atmosphere, but also on its temperature structure. The influence of the decreasing stellar 
effective temperature of M0-M7 stars on the planetary atmos pheres can be seen in our simulated 
planet spectra. As already pointed out bv lSegura et al.l (12005b . absorption of O3, CH4, and N2O is 
easier to detect for planets that orbit active M dwarfs. We confirm this result for AD Leo and also 
for quiet M-dwarf central stars up to about type M3. For cooler central stars, however, the develop- 
ing temperature inversion in the mid-atmosphere leads to reduced absorption signals, which almost 
level out for the coolest star considered (M7). Thus, although the abundances for most molecules 
continuously increase from M0 to M7 central stars, absorption features of atmospheric molecules 
are very difficult to detect in emission spectroscopy of planets that orbit very cool M dwarfs, which 
lead to potentially false negative detections of biomarkers in these cases. 

Our simulations of transmission spectra during primary eclipses show that the relative trans- 
mission in spectral absorption bands of CO2, H2O, CH4 increases from M0 to M7 dwarfs, hence 
increases towards cooler M central stars. We point out that the increasing relative transmission of 
the biomarker N2O becomes comparable with the major CO2 band for Earth-sized planets around 
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M7 stars. Because N2O is mainly produced from denitrifying bacteria on Earth, it is considered to 
be a good biomarker molecule. Unfortunately, for most planet scenarios the N2O absorption bands 
are far too weak to be detected compared with the dominant atmospheric constituents. This may be 
different for planets around central stars with very low UV fluxes . 

The situation is different for transmission spectra of the biomarker ozone. For active and hot 
M dwarfs ozone absorption also increases. However, for planets around cool M stars (M4-M7) the 
ozone absorption is reduced again because of its smaller atmospheric abundances caused by the 
low UV flux in quiet M dwarf stars. 

Our simulations confirm that low spectral resolution (around R=20) is sufficient to identify 
most spectral features in the spectrum of a planet with an Earth-like atmosphere for IR emission 
as well as transmission spectra. A characterisation of an atmospheric temperature structure using 
IR emission spectra, e.g. by studying the presence of a temperature inversion, would require much 
higher spectral resolutions (of about >200). 

Generally, emission spectra are studied in the IR wavelength range. Our SNR calculations show 
that the detection of the broad CO2 absorption feature at 15 //m and the ozone band at 9.6 fim could 
be feasible for planets around very nearby M dwarfs with exposure times of several hours for our 
modelled planet scenarios. 

Transmission spectra obtained during primary eclipse of the planet, on the other hand, favour 
shorter wavelength ranges owing to the increasing stellar flux and therefore better SNR. In this 
case, detection of e.g. methane at 3.3 fim may be possible when adding transits for the nearest 
M dwarfs. Summing many transits with a stable instrument could bring the measurements into a 
feasible range for future observations. 

Of course, the real atmospheric composition and total pressures of super-Earth planets remain 
unknown until the first spectroscopic detections can be made. Atmospheres even of planets sim- 
ilar to Earth may differ in atmospheric mixing ratios, or show e.g. more extended atmospheres . 
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2010). 



Clouds will additionally affect the detection of absorption signals (see e.g 
This will result in different SNR from those computed here, where only particular scenarios were 
studied. Investigations of a wide parameter range would therefore be needed. Looking for yet un- 
detected planets, it is important to be open minded and take into account the unexpected, e.g. 
when designing instruments. Consistent modelling of climate and chemistry in the atmospheres of 
extra-solar planets is important, even when possible only in exemplary cases, for understanding 
effects on T-p-profiles and on the composition and therefore spectral appearance. An example is 
the effect of the T-p-profiles for emission spectra of cool M dwarf stars, which in the extreme case 
leads to extremely shallow absorption signals and consequently false negative detections. Another 
example is the ozone signal in transmission spectra, where the ozone column is significantly re- 
duced for planets around very cool stars owing to chemical effects. For the main CO2 15 fim band, 
however, the differences are small when only considering detection of the broad absorption band. 
Nevertheless, if future telescopes will be able to observe terrestrial exoplanets at higher spectral 
resolution (R>100), a consistent model approach for CO2 is also mandatory to interpret the tem- 
perature inversion peak of this band. 
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